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Abstract - The consequences of dynamical screening of Coulomb interaction among correlated 
electrons in realistic materials have not been widely considered before. In this letter we try to in- 
corporate a frequency dependent Coulomb interaction into the state-of-the-art ab initio electronic 
structure computing framework of local density approximation plus dynamical mean-field theory, 
and then choose SrVC>3 as a prototype material to demonstrate the importance of dynamical 
screening effect. It is shown to renormalise the spectral weight near the Fermi level, to increase 
the effective mass, and to suppress the t^ g quasiparticle band width apparently. The calculated 
results are in accordance with very recent angle-resolved photoemission spectroscopy experiments 
and Bose factor ansatz calculations. 



Introduction. — Over the last years, tremendous pro- 
gresses have been made in the investigations of strongly 
correlated materials. [IH3] Such systems typically contain 
partially filled and relatively localized 3c?, 4/ or 5/ or- 
bitals, and the low-energy physics of these systems is gen- 
erally described by an effective Hamiltonian which mod- 
els the electronic Coulomb interaction among these cor- 
related orbitals. The electronic Coulomb interaction can 
be generally parameterized by a number, the "Hubbard 
U" . Usually U is assumed to be static, but in principle 
U is dynamical and frequency dependent: electric fields 
will be induced by the charge density fluctuation in cor- 
related orbitals, while the higher or lower-lying states in 
the systems will act to screen, resulting in a frequency 
dependent renormalization. [1HS] On one hand, it is real- 
ized that at high frequency the screening becomes weaker 
and eventually the interaction approaches the large bare 
Coulomb value, which is an order of magnitude larger than 
the static screened value. On the other hand, it is gener- 
ally believed that the dynamical screening effect may play 
a key role in understanding the subtle electronic struc- 
tures of strongly correlated materials, [?H9] for examples, 
plasmon satellites and spectral weight transfer in photoe- 
mission spectra of hole-doped BaFe2As2. jTOj Thus how to 
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take the dynamical screening effect into consideration is 
one of the challenges in studying the intriguing properties 
of strongly correlated materials. 

In a realistic approach to strongly correlated materials, 
the frequency dependent U = U (u) can be evaluated at 
the random phase approximation (RPA) level, which al- 
lows for the direct evaluations of matrix elements of Hub- 
bard U and its frequency dependence. [IHS] The static 
screened value Uq = 5RJ7(o; = 0), has been widely used 
in the local density approximation plus static Coulomb 
interaction (LDA + U) method [TTJ[T2] and local den- 
sity approximation combined with dynamical mean-field 
theory (LDA + DMFT) approach. [^[IMS] Since the 
frequency dependent screening has been completely ne- 
glected or empirically taken into account by adjusting the 
effective static Uq, very little is known on the consequences 
of frequency dependent Coulomb interaction in the elec- 
tronic structures of strongly correlated materials. In order 
to include the dynamical screening effect into the LDA 
+ DMFT framework, the hardest obstacle has been the 
lack of a reliable and efficient quantum impurity solver for 
the general impurity model with a frequency dependent 
Coulomb interaction U(u>). [5] 

Fortunately, this obstacle has been recently overcome 
by the development of strong coupling continuous-time 
quantum Monte Carlo impurity solver (abbreviated CT- 
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HYB) proposed by Werner et al. P ITMTB] . where a multi- 
plasmon Lang-Firsov transformation is treated exactly in 
the context of a hybridization expansion algorithm for the 
general impurity model. By using this powerful impurity 
solver, Werner et al. have studied the plasmon satellites 
and large doping and temperature dependence of elec- 
tronic properties in hole-doped BaFc2As2, a fascinating 
iron-based superconductor. They have demonstrated that 
the dynamical screening effect is important not only for 
high-energy spectral features, such as correlation satel- 
lites seen in photoemission spectroscopy, but also for the 
low-energy electronic structure, mass enhancements and 
lifetimes of hole-doped BaFe2As2- [10] As an very useful 
complement, the weak coupling continuous-time quantum 
Monte Carlo impurity solver (abbreviated CT-INT) devel- 
oped by Rubtsov et al. [TM2T)] can treat generic retarded 
interactions as well. But the CT-INT impurity solver 
is limited to small dynamical screened U(u) and not-so- 
large screening frequencies ui, |20j and therefore becomes 
prohibitively costly for realistic applications of dynamical 
screening interactions in multi-orbital impurity models. 
M. Casula et al. [9] proposed a dubbed Bose factor ansatz 
(BFA) for the Green's function of quantum impurity prob- 
lem with retarded interactions, in which the Green's func- 
tion is factorized into a contribution stemming from an 
effective static-C/ problem and a bosonic high-energy part 
introducing collective plasmon excitations. Various ap- 
proximations for the Green's function Bose factor F are 
introduced and their pros and cons arc analyzed care- 
fully. The most practical and effective one is borrowed 
directly from the so-called dynamical atomic limit ap- 
proximation (DAL A), whose functional form is analyti- 
cally known. Since the DALA is valid in the hybridiza- 
tion A = limit, it only works well in the intermediate- 
strong coupling regime, with Uq and the dynamical part 
large enough. The other functional forms going beyond 
the DALA scheme are available as well. [9] 

Inspired by Werner et al.'s pioneer works, OHO] in this 
letter we try to incorporate the dynamical screening effect 
into the LDA + DMFT computing framework by adopt- 
ing their approach. [8]IT6UT8] Then we apply this technique 
to solve a realistic multi-band Hamiltonian with dynami- 
cal Coulomb interaction U(uj) for prototypical transition 
metal oxide SrV03. The quasiparticle weight reduction, 
more renormalized tig band width and the spectral weight 
transfer to higher frequencies are found and analyzed in 
details. Our calculations highlight the importance of in- 
cluding the proper screened interaction, to have a cor- 
rect and reliable description of the electronic correlation 
in strongly correlated materials. 

Method. — The computational scheme used in 
present work can be viewed as an extension of the tra- 
ditional LDA + DMFT method. [ffl3] In a Hamiltonian 
formulation, the multi-orbital model with dynamically 
screened interactions can be written as 



Here -ff™ A is the single particle low-energy effective 
Hamiltonian. To generate -f^DA, the ground state cal- 
culation by using the pseudopotential plane- wave method 
with the abinit package [21] has been performed at first. 
For the ground state calculation, the projector augmented 
wave (PAW) type pseudopotentials [22] for Sr, V and O 
species are constructed by ourselves, the cutoff energy for 
plane-wave expansion is 20 Ha, and the fc-mcsh for Bril- 
louin zone integration isl2xl2xl2. These pseudopoten- 
tials and computational parameters are carefully checked 
and tuned to ensure the numerical convergences. And 
then -f^LDA i s obtained by applying a projection onto max- 
imally localized Wannier function (MLWF) [23] orbitals 
including all the Y-t 2 g and 0-2p orbitals, which result- 
ing in a 12 x 12 tight-binding Hamiltonian. H-qq is the 
double counting term, and the around mean-field (AMF) 
scheme [T5] is used, which is especially suitable for strongly 
correlated metal system. The Coulomb interaction term 
Hi is taken into consideration merely among three V-i2g 
orbitals. In the present work, we choose U = 4.0 and J = 
0.65 eV, which are close to previous estimations. [13, 15 
The bosonic part Hs reads 



Ho 



, (2) 



which describes the coupling of electronic degrees of free- 
dom (occupation number n) to bosonic modes (bosonic 
operators and 6) and the coupling strength \ u = 

^sll (u>) 1 7T . The bosonic modes represent the screening of 
the charge density fluctuations on the strongly correlated 
V-^2g orbitals. 

The multi-orbital lattice problem is studied using 
DMFT, [TJE] which maps it to a self-consistent solution of 
a three-orbital quantum impurity model. Then the impu- 
rity model (with frequency-dependent interactions U(u)) 
is solved using the CT-HYB impurity solver 8, 16 18 
which based on a stochastic diagrammatic expansion of 
the partition function in the impurity-bath hybridization. 
In fact, the effect of U{uj) is to dress the fermionic propa- 
gators with a bosonic propagator cxp [— /C(r)] where /C(t) 
is the twice-integrated retarded interaction. In terms of 
$sU(u>) and a factor B(t,lj) = cosh[(r — §)w]/ sinh[4p] 
with bosonic symmetry, /C(r) can be expressed by [8]ll0j 



/C(r) = 



dui 5sU(ui) 



[B(t,u) - B(0,u>)]. 



(3) 



H 



Hi + H s . 



(1) 



Thus the diagrammatic weight for Monte Carlo stochas- 
tic sampling is supplemented with an additional bosonic 
term, which originates from the contribution of /C(r). The 
implemented details of CT-HYB impurity solver improved 
with recently developed orthogonal polynomial represen- 
tation algorithm [24] can be easily found in the litera- 
tures. [SQI1QI] In each LDA + DMFT iterations, typi- 
cally 4 x 10 8 Monte Carlo samplings have been performed 
to reach sufficient numerical accuracy. All the calcula- 
tions are carried out at the inverse temperature of /3 = 10 
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Fig. 1: (Color online) Frequency dependent interaction for 
SrVC>3 from constrained random phase approximation, (a) 
Real and imaginary parts of the average intra-orbital Coulomb 
interaction as a function of frequency. The data shown in this 
figure are taken from reference [5j directly, (b) The twice- 
integrated retarded interaction function /C(r) at finite temper- 
atures j3 = 10 and /3 = 30 calculated by Eq.([3]). 



or /3 = 30. Finally the maximum entropy method |25j 
is used to perform analytical continuation to obtain the 
impurity spectral functions of V-^g states. 

Results and discussion. — SLVO3 is a prototype 
of correlated metal and very well studied material which 
represents a benchmark for theories describing strongly 
correlated compounds. Indeed, its band structure is 
relatively simple due to its undistortcd pcrovskite struc- 
ture, resulting in the occupation of one electron in three- 
fold degenerate W-t 2g bands crossing the Fermi level. The 
Oxygen 2p bands are quite well separated from the t 2g 
levels, such that the 3x3 V-i2 g tight-binding Hamilto- 
nian is a minimal model required for a correct description 
of the low-energy physics of SrVC>3. Thus, SrVC>3 has 
been the testing case for many new LDA + DMFT imple- 
mentations. [TMUIMSn] And on the other hand, SrV0 3 
has been the subject of intensive experimental activ- 
ity, with magnetic, electrical, optical measurements, and 
by means of angle-resolved photocmission spectroscopy 
(ARPES). [3DHS3] Therefore SrV0 3 is an ideal system to 
benchmark our newly developed implementation of LDA 
+ DMFT computing framework incorporated with dy- 
namical screening interaction. 

To include dynamical screening interaction within LDA 
+ DMFT computing scheme, one of the bottlenecks is the 
determination of the frequency dependent Coulomb inter- 
action U (w) and corresponding twice-integrated retarded 
interaction /C(t). Recent works try to extract these quan- 
tities from constrained density functional calculations, [TT] 
GW-inspircd methods, [3] or RPA-based schemes. [SJ|S] 
All these calculations are very tedious and tricky. Fortu- 
nately, the realistic frequency dependent U(uj) for SrV03 
has been computed before [5] based on the constrained 
RPA (cRPA) approach, so we can extract U(ui) directly 
from the reference as a necessary input to evaluate JC(t) 
in the preprocessing stage. The obtained U(oj) and /C(r) 
are shown in FigfT] respectively. Clearly seen in FigQTa), 



static U 

dynamic U 

(a) ^^J)\ 








IN^-^^^Vv -20 -15 -10 -5 15 20 25 30 


(b) ^_Jr 




A 


A 


-20 -15 -10 -5 15 20 25 30 


(c) _J 








" ,s: ^ > ^ -20 -15 -10 -5 15 20 25 30 



5 10 
Frequency (eV) 



0.45 
0.4 

0.35 

0.3 

& 0.25 
- 

0.15 
0.1 
0.05 










c 


DALA -e- 

T-HYB - 






















< 



























































2 
Frequency (eV) 



Fig. 2: (Color online) Integrated spectral functions A(uj) of 
V-t2 9 states. Upper panel: Spectral functions obtained in 
current LDA + DMFT calculations for different models, (a) 
iVwann = 3, j3 = 10; (b) iVwann = 3, p = 30; (c) iV wa „„ = 
12, /3 — 10. Inset: The fine structures of spectral func- 
tions at high frequencies. The spectral functions are obtained 
from imaginary-time Green's functions G(t) by using maxi- 
mum entropy method, [55] and the results are cross-checked 
by using recently developed stochastic analytical continuation 
method. [34] Lower panel: Comparison of spectral functions of 
V-t2 9 states obtained by our LDA + DMFT calculations with 
CT-HYB impurity solver and BFA-DALA method respectively. 
iVwann = 3, @ = 10. The results obtained by BFA-DALA 
method are extracted directly from the reference [5]. 



^sU(uj) features prominent peaks at 14 eV, 26 eV and 33 
eV respectively. From the peaks of SC/(w)/u; 2 , the ener- 
gies of plasmon satellite structures in the V-^2g density 
of states can be easily concluded. The real part of U(oj) 
ranges from the static screened value Uq = 4.0 eV to the 
bare unscreened value of about 12 eV at large oj, and the 
frequency dependence resembles typical U(uj) in strongly 
correlated metals. In a standard LDA + DMFT calcu- 
lation without dynamical screening effect, the relatively 
small value of static screened value Uq would result in a 
rather weakly correlated picture. 

Now let's focus on the spectral properties of V-i2g 
states. Two different models for SrV03 are taken into 
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considerations, one contains only the V-i2g states (num- 
ber of Wannier orbitals iV wann = 3), the other contains 
both V-t2 g and 0-2p states (number of Wannier orbitals 
JV WMm = 12). We solve these models in the framework 
of LDA + DMFT at = 10 and = 30 with dynam- 
ical U and static U respectively. The calculated results 
are illustrated in the upper panel of FigJ5J The results 
obtained by dynamical U and static U calculations dis- 
play remarkable differences. Firstly in the dynamical U 
calculations, spectral weights arc shifted strongly to high 
frequencies, which leads to a reduction in weight at low 
energies, compared with the static U results. Secondly, 
the quasiparticle resonance peak around the Fermi level 
is strongly renormalized, not only the width but also the 
height of it shrink apparently as U(uj) is taken into ac- 
count. Thirdly, in dynamical U calculations, some shoul- 
der peaks near w = arc smeared out. Thus, the explicit 
treatment of the Coulomb repulsion at large frequencies 
has a substantial effect, even on the low-energy properties 
of the system. [MtO] 

We note that recently Casula et al. [9] have studied the 
spectral properties of SrVC>3 by using their BFA-DALA 
method. Since both the BFA-DALA method and our CT- 
HYB scheme are intent to deal with the strongly correlated 
systems with frequency dependent Coulomb interaction, it 
should be meaningful to justify the calculated results ob- 
tained by these two different methods. Comparison of the 
calculated spectral functions of V-t i g states arc 1 shown in 
the lower panel of Figj2] Both spectral functions exhibit 
significant three-peak structures, which are prominent for 
strongly correlated metal systems, p] Clearly, the spec- 
tral function obtained by our CT-HYB scheme is consis- 
tent with that obtained by BFA-DALA method, besides a 
slightly spectral weight reduction near the Fermi level and 
a small shift of the upper Hubbard band to high energies. 
These discrepancies are certainly due to the smaller value 
of static screened U (Uq = 3.6 eV) used in Casula et al.'s 
BFA-DALA calculations. [9] 

In Figure [3] the calculated self-energy function of V-t2 g 
states are reported. Let's concentrate our attention to the 
imaginary part of self-energy function at matsubara fre- 
quency axis at first. In dynamical U model, — Sy£(iu; — »• 0) 
is remarkable larger than that of corresponding static U 
model. According to the well-known Eliashberg equa- 
tion: [1] Z^ 1 wl- £$sH(iuj — > 0), the quasiparticle weight 
Z can be approximately evaluated. Thus the dynamical 
screening effect results in a smaller value of Z, namely a 
larger value of effective mass m* . Indeed, for the realistic 
dynamical U model we found a Z w 0.43 ~ 0.48, which 
gives an effective mass renormalized by 2.1 ~ 2.3 with re- 
spect to the LDA band structure. On the other hand, for 
the corresponding static U model we obtained a value of 
Z w 0.6, which underestimates the electronic correlation 
by a factor of 1.67, and so the value of m* . Recent ARPES 
data yielded an effective mass m* « 2mo, [3"0H3"2"] which is 
in good agreement with our findings for the realistic dy- 
namical interaction. It is worthy noting that a static U 
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Fig. 3: (Color online) Self-energy function of V-t2d states ob- 
tained by LDA + DMFT calculations at finite temperature. 
Left paneklmaginary part of self-energy function at matsubara 
frequency axis QY*(iuj). Right panelReal part of self-energy 
function at real axis KE(oj). The blue and green arrows are 
used to indicated the "kink"-like structures, which are often 
considered as signatures of strongly correlated systems. 



model with a larger instantaneous Uq [26H28] can be used 
to fit the experimental mass renormalization artificially. 
But the lower Hubbard band peaked at —1.5 eV can not 
be obtained by such a static model. The difficulty here 
is to reproduce both the effective mass and the position 
of the lower Hubbard band by the same model. With the 
dynamically screened interaction model, we can describe 
correctly both the effective mass (~ 2.1mo) and the lower 
Hubbard band peaked at ~ —1.5 eV, as is seen in both 
Fig! and FigEI 

After careful analytical continuation, self-energy func- 
tion of V-i2 S states at real axis is obtained and plot- 
ted in Fig[3]Jb). For the real part of self-energy function 
3?S(w), there are two extrema at the energies u> ~ ±1.0 
eV, originating from the crossover from the central quasi- 
particle resonance peak to the lower and upper Hubbard 
bands. In the energy regime of the quasiparticle reso- 
nance peak ranging from about -0.5 eV to 0.5 eV, 3£E(uj) 
can be roughly described by a straight line and the slope 
of it can be used to evaluate the quasiparticle weight Z as 
well. Not surprisingly, the quasiparticle weight Z calcu- 
lated with 3ftS(w) is almost identical with that calculated 
with S5S(zu;). In this energy range the Landau Fermi liq- 
uid theory is valid. Strictly speaking, the Fermi liquid 
regime with KE(cj) ~ — tu only extends from -0.2 up to 
0.2 eV, which is in accord with previous LDA + DMFT 
calculations. [26] Next to this Fermi liquid regime, there 
are pronounced "kink"-like structures in KS(w) around 
uj = ±0.25 eV as are highlighted by colourful arrows. Ap- 
parently, the "kink" -like structures are enhanced by the 
dynamical screening effects. These "kink" -like structures 
can be regarded as a fingerprint of strongly correlated sys- 
tems |35) and become important in the context of quasi- 
particle dispersion and electronic specific heat. In previ- 
ous LDA + DMFT calculation and ARPES experi- 
ment [30] the "kink" -like structures in SrV03 have been 
already confirmed. 
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Fig. 4: (Color online) Angle-resolved spectral functions A(k, lj) 
for SrVC>3. Upper panel: Static Coulomb interaction case. 
Lower panel: Frequency dependent Coulomb interaction case. 
The colourful dash lines denote LDA band structures and the 
green arrows indicate "kink"-like structures. In these LDA 
+ DMFT calculations, only three V-t2 a bands are taken into 
considerations and the inverse temperature /3 is fixed to 10. 



With the knowledge of self-energy function on the real 
axis, we are now in the position to calculate k-resolvcd 
spectral functions or quasiparticle dispersion A(k,uj). The 
LDA + DMFT quasiparticle dispersions for dynamical U 
and static U models along high symmetry lines in the Bril- 
louin zone are shown in Fig|2] respectively. The LDA band 
structure is presented in this figure as an useful compar- 
ison. As for the static U model, the quasiparticle band 
structure displays strongly band renormalization with a 
factor ~ 1.7 with respect to LDA band dispersion. For 
example, we can see from the upper panel of FigQ] that 
the bottom of the quasiparticle band is located at approx- 
imately u) = —0.6 cV, in contrast to the LDA value of 
w = —1.0 cV. Near the Fermi level along R — T and M — T 
lines there exist discernible "kink" -like structures as in- 
dicated by a green arrow, [551[3D] which stem from the 
shoulders in the real part of the self-energy function (see 
FigCUb)). As for the dynamical U model, its quasiparti- 
cle band structure exhibits similar features with respect 
to that of static U model. Around the Fermi level, the 



quasiparticle band structure of dynamical U model is more 
renormalized with a factor ~ 2.1 and has less intensity, 
and the bottom of quasiparticle dispersion is located at 
lu = —0.5 eV. The "kink" structures in quasiparticle band 
structure are more discernible than those of corresponding 
static U model. 

Conclusion. — In this letter, we try to incorporate the 
dynamical screening effect into the modern LDA + DMFT 
computing framework with CT-HYB as a quantum impu- 
rity solver. We apply this calculation scheme to reinves- 
tigate the electronic structure of SrV03 with dynamical 
U model. It seems that the dynamical screening effect 
will enhance the electron correlation significantly. The 
transfer of spectral weight to higher frequencies, further 
reduction of quasiparticle resonance peak, smaller quais- 
particle weight Z, larger effective mass m*, more apparent 
"kink" -like structures etc. are predicted or verified by our 
calculations as well. Based on the calculated results, it 
is suggested that the dynamical screening interaction may 
play an important role in understanding the fine electronic 
structures of strongly correlated materials, which has been 
ignored by most of previous theoretical calculations. 
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